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Abstract. By studying the present cosmological data, particularly on CMB, SNelA and LSS, we 
find that the future fate of the universe, for simple linear models of the dark energy equation-of- 
state, can vary between the extremes of (I) a divergence of the scale factor in as little as 7 Gyr; (II) 
an infinite lifetime of the universe with dark energy dominant for all future time; (III) a disappearing 
dark energy where the universe asymptotes as f ^ oo to a(f ) ^ f^/^ i.e. matter domination. Precision 
cosmological data hint that a dark energy with equation of state w — P/p < — 1 and hence dubious 
stability is viable. Here we discuss for any w nucleation from A>0toA = 0ina first-order phase 
transition. The critical radius is argued to be at least of galactic size and the corresponding nucleation 
rate glacial, thus underwriting the dark energy's stability and rendering remote any microscopic 
effect. 



INTRODUCTION 



In this talk I will focus on two issues concerning the enigmatic dark energy which 
comprises approximately 70% of the cosmological energy density in the present era. 

The first is concerning the future fate of the dark energy and how it depends on the 
equation of state. 

The second is the issue of stability of the dark energy assuming that it can undergo a 
first-order phase transition to a vacuum of truly zero energy density. 



FATE OF DARK ENERGY 



The cosmic concordance of data from three disparate sources: Cosmic Microwave Back- 
ground (CMB), Large Scale Structure (LSS) and High-Red-Shift Supernovae (SNelA) 
suggests that the present values of the dark energy and matter components, in terms of 
the critical density, are approximately D.x — 0.7 and Q.m — 0.3. The question to which 
we try to make a small contribution in this paper is to what extent precision cosmological 
data will allow us to discriminate between possible futute fates of the Universe? 

If one makes the most conservative assumption that Q.x corresponds to a cosmological 
constant with Equation of State given by a constant w = p/p = —I, then the future 
evolution of the universe follows from the Friedmann equation 
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in which a{t) is the scale factor normalized at the present time as a{to) — I, p{t) = 
p (0)a~^ is the energy density of matter component and A is constant. Here we assumed 
that the universe is flat and neglected radiation. 

In such a simple, and still viable, case the behavior of a{t) for asymptotically large 

r — > oo is 



so that the dark energy asymptotically dominates and the universe is blown apart in an 
infinite time a{t) — > <» as ? — > <». 

Even assuming that w is constant, however, there is a wide range of possible w: 
according to [1] the allowed values are —2.68 < w < —0.78. We do not assume this 
result but will arrive at a similar allowed range; the difference is because our priors are 
slightly different (we fix the cosmological parameters CIm = 0.3, Q.\ = 0.1, Q.^ = 0.02, 
and h = 0.65 instead of allowing them to vary.) The asymptotic behavior in Eq.(2) 
is very far from established by present data. The case w < — 1 has the property that 
boosting from the dark energy rest frame to an inertial frame with velocity satisfying 
(v/c)^ > — 1/w leads to a negative energy density, but this does not violate any law 
of physics^ Here we shall consider some simple models for w, including dependence 
on red-shift w{z}, to illustrate how far existing data are from answering the question 
of the future fate of the universe. As we will show, for a model in which w varies 
linearly with red-shift, present data are consistent with extremely different futures. For 
examples, in one case the scale factor diverges[4, 5] in finite time ^, in just another 7 
Gyr, while in another case the energy density of dark energy decreases eventually faster 
than that of matter, i.e., the dark energy disappears and the universe reverts to being 
matter-dominated with a{t) ~ t^^^, as ? — > <». 

Constant Equation of State. 

Here we discuss the future fate of the universe in the case of the constant equation of 
state. If we assume, to begin, that w is constant then, keeping only the dark energy term 



where j8 = 3(1 Most authors have discussed the case with w > — 1, however the 
case with w < — 1 is also possible and discussed phenomenologically in [1, 4, 6], and in 



This violates the weak energy condition[2, 3]. 
^ Gravitationally-bound systems could survive longer than tr in Eq.(5) but such systems would be in- 
finitely separated from one another. 
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connection with string theory in [7]. If /3 < 0, corresponding to w < —1, the solution of 
Eq.(3) diverges at a finite time t — t*.By integrating 

poo pt 

/ a^l^-'^Ho^ / dt (4) 

one finds that the remaining time tr before time ends tr = {t* — to) is given analytically 
by 
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In Eq.(5), putting in Qx = 0.1 and ^H^^ = 9.2 Gyr one finds for w = -1.5, -2.0 and 
—2.5, respectively ti = 22, 11 and 7.3 Gyr. 

In such a constant w scenario which is consistent with all cosmological data, the 
divergence of the scale factor will occur in a finite time period of 7 Gyr (or more) from 
now. 

With respect to the Solar System, this end of time occurs generally after the Sun has 
transformed into a Red Giant, and swallowed the Earth, as is expected approximately 5 
Gyr in the future. 



Equation of State Varying Linearly with Red-Shift. 

To view the four Figures cited, see [25] 

As a more general ansatz, we consider the model for w depending linearly on red- 
shift:^ 

w(z) = w(o)+cze{C-z)+cce{z-Q (6) 

where the modification is cut off arbitrarily at some Z — ^ > 0. We assume C < and 
consider the two-dimensional parameter space spanned by the two variables w(0) and 
C. 

In order to discuss constraints on our phenomenological model, we compare its 
prediction with experimental data from SNelA and CMB. We evaluate the goodness-of- 
fit parameter x'^ as a function of w{0) and C. For SNelA, we used a dataset consisting 
of 37 SNe from [10] with the MLCS method. To calculate the CMB power spectrum, 
we used a modified version of CMBFAST [11]. For the analysis of CMB, we used 
experimental data from COBE [12], BOOMERanG [13], MAXIMA [14] and DASI [15]. 
To calculate x'^, we adopt the offset log-normal approximation [16] and used RADPACK 
package [17]. We also studied the constraint from LSS using the 2dF data [18], and found 
that it does not give severe constraint on the parameters w{0) and C. 



^ A model with w linearly depending on red-shift is also discussed in [8] but fitting the CMB data was not 
investigated. Another parametrization of w(Z) is in [9]. 



For an illustration, we take the cosmological parameters as Q.\)h^ = 0.02, Qm = 
0.3, = 0-7 and h — 0.65, and the initial power spectra are assumed to be scale 
invariant in all numerical calculation in this paper. 

To set the stage, let us first use only the SNelA data to constrain the parameters w(0) 
and C. The result is shown for ^ = 2 in Figure 1 where the 99 % C. L. allowed region is 
the region between the two dashed lines shown. We may remark three distinct regions: 

(I) w{0) < (C — 1). In this case there is divergence of the scale factor, at a finite future 
time. 

(II) (C — 1) < w{0) < C. Here the lifetime of the universe is infinite. The dark energy 
dominates over matter, as now, at all future times. 

(III) C <w{0). The lifetime of the universe is again infinite but after a finite time the 
dark energy will disappear relative to the dark matter and matter-domination will be 
re-established with ~ 

When we add the constraints imposed by the CMB data, the allowed region is smaller 
as shown in Figure 2, plotted for ^ = 0.5. Such a small ^ still allows all three future 
possibilities (I), (II) and (HI). For somewhat larger ^ only possibilities (I) and (II) are 
allowed in this particular parameterization. 

The case ^ = 2 is exhibited in more detail for different values of h'(O) and C in Figures 
3 and 4. Figure 3 shows the variation of the transition red-shift Ztr where deceleration 
changes to accelerated cosmic expansion defined by q{Z^) = 0. From the figure, we can 
read off that Ztr becomes smaller as C becomes more negative for fixed w(0); this is 
because the epoch where dark energy becomes the dominant component of the Universe 
becomes later. This affects the magnitude-red shift relation of high-Z supemovae. In 
Figure 4, the magnitude-red shift relation is shown for the SNelA data [10] along with 
the prediction of our phenomenological model for ^ = 2. The magnitudes are calculated, 
as usual, relative to the empty universe Milne model with D.m — 0, D.x — and Q.j^~ I. 
One can expect that high-redshift SNe would appear dimmer if C were more negative. 

To return to our main point, let us assume that more precise cosmological data will 
allow an approximate determination of w{Z) = f{Z) as a function of Z for positive Z > 0. 
Then to illustrate the possible future evolutions write: 

w(Z) = /(Z)0(Z) + (/(O) + aZ)0(-Z) (7) 

In this case, the future scenarios (I), (II) and (IE) occur respectively for a > (/(O) + 1), 
(/(0) + l)>a>/(0)anda</(0). 

Present data are consistent with a simple cosmological constant /(Z) = — 1 in Eq.(7) 
in which case the divergence of the scale factor occurs for a > 0, the infinite-time dark 
energy domination for > a > —1, and disappearing dark energy for a < —1. 

Since in practice F{Z) for Z > will never be determined with perfect accuracy 
the continuation of w(Z) to future Z < will be undecidable from observation as will 
therefore be the ultimate fate of the Universe. 



STABILITY ISSUES FOR DARK ENERGY 



The equation of state for the dark energy component in cosmology has been the subject 
of much recent discussion [7, 19, 20, 21, 22, 23, 24, 25, 26] Present data are consistent 
with a constant w{Z) = — 1 corresponding to a cosmological constant. But the data allow 
a present value for w{Z = 0) in the range —1.38 < w{Z = 0) < —0.82 [19]. If one 
assumes, more generally, that w(Z) depends on Z then the allowed range for w{Z — 0) is 
approximately the same[25]. In the present article we shall forgo this greater generality 
as not relevant. Instead, in the present article we address the question of stability for a 
dark energy with constant w(Z) < —1. 



Interpretation as a limiting velocity 



Consider making a Lorentz boost along the 1 -direction with velocity V (put c : 
the stress-energy tensor which in the dark energy rest frame has the form: 



r^v = A 
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is boosted to T^^ given by 
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We learn several things by studying Eq.(9). First, consider the energy component Tqq = 
1 -|- V w. Since V < 1 we see that for w > — 1 this is positive Tqq > and the Weak 
Energy Condition (WEC) is respected[27]. For w = — 1, Tqq — > as V ^ 1 and is still 
never negative. For w < — 1, however, we see that Tqq < if y2 > and this 

violates the WEC and is the first sign that the case w < —I must be studied with great 
care. Looking at the pressure component T^^ we see the special role of the case w = — I 
because w=T^i/Tqq remains Lorentz invariant as expected for a cosmological constant. 
Similarly the off-diagonal components Tqj remain vanishing only in this case. The main 
concern is the negativity of Tqq < which appears for > — (1/w). One possibility 
is that it is impossible for > — (1/w). The highest velocities known are those for 
the highest-energy cosmic rays which are protons with energy ~ 10-^^ eV. These have 
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which is one possible conclusion. 



First-Order Phase Transition and Nucleation Rate. 

But let us suppose, as hinted at by [19, 20, 21] that more precise cosmological data 
reveals a dark matter which violates Eq.(lO). Then, by boosting to an inertial frame 
with > — (1/w), one arrives at Tqq < and this would be a signal for vacuum 
instability [30]. If the cosmological background is a Friedmann-Robertson- Walker 
(FRW) metric the physics is Lorentz invariant and so one should be able to see evidence 
for the instability already in the preferred frame where T^v is given by Eq.(8). This goes 
back to work in the 1960's and 1970's where one compares the unstable vacuum to a 
superheated liquid. As an example, at one atmospheric pressure water can be heated 
carefully to above 100^ C without boiling. The superheated water is metastable and 
attempts to nucleate bubbles containing steam. However, there is an energy balance 
for a three-dimensional bubble between the positive surface energy ~ and the 
negative latent heat energy of the interior ~ R^ which leads to a critical radius below 
which the bubble shrinks away and above which the bubble expands and precipitates 
boiling[28, 29]. For the vacuum the first idea in [30] was to treat the spacetime vacuum 
as a four-dimensional material medium just like superheated water. The second idea in 
the same paper was to notice that a hyperspherical bubble expanding at the speed of light 
is the same to all inertial observers. This Lorentz invariance provided the mathematical 
relationship between the lifetime for unstable vacuum decay and the critical radius of 
the four-dimensional bubble or instanton. 
In the rest frame, the energy density is 

roo = A- (10"^ey)^~ (Imm)^"^ (11) 

since lO^^eV ~ (Imm)^^. 

In order to make an estimate of the dark energy decay lifetime in the absence of a 
known potential, we can proceed by assuming (without motivation from observation) 
that there is a first-order phase transition possible between the A — (10~^eV')'* "phase" 
and a A = "phase" . This hypothesized decay is the Lorentz invariant process of a 
hyperspherical bubble expanding at the speed of light, the same for all inertial observers. 
Let the radius of this hypersphere be R, its energy density be e and its surface tension 
be ^1. Then according to [30] the relevant instanton action is 

A = -]^7t^R'^e + lTi^R^Si (12) 

where e and are the volume and surface energy densities, respectively. The stationary 
value of this action is 

Am=^-Tt'S!\/e^ (13) 

corresponding to the critical radius 



Rm = 3Si/e 



(14) 



We shall assume that the wall thickness is negligible compared to the bubble radius. The 
number of vacuum nucleations in the past lightcone is estimated as 

N = {VuA'^)expi-A,n) (15) 

where Vu is the 4-volume of the past and A is the mass scale relevant to the problem. 
This vacuum decay picture led to the proposals of inflation[31], for solving the horizon, 
flatness and monopole problems (only the horizon problem was generally known at 
the time of [30]). None of that work addressed why the true vacuum has zero energy. 
Now that the observed vacuum has non-zero energy density +e ~ {lO~^eV)'^ we may 
interpret it as a false vacuum lying above the true vacuum with e = 0. In order to use the 
full power of Eq.(15) taken from [30, 32] and the requirement A'^ ^ 1 we need to estimate 

the three mass-dimension parameters £^/'^,s\^^ and A therein and so we discuss these 
three scales in turn. 

The easiest of the three to select is e. If we imagine a tunneling through a barrier 
between a false vacuum with energy density £ to a true vacuum at energy density zero 
then the energy density inside the bubble will be £ = A = (lO^^eV)^. No other choice 
is reasonable. As for the typical mass scale A in the prefactor of Eq. (15), the value of 
A does not matter very much because it appears in a power rather than an exponential 
so let us put (the reader can check that the conclusions do not depend on this choice) 
A = £^/^ = (Imm)^^ whereupon the prefactor in Eq.(15) is ~ 10^^^. The third and final 
scale to discuss is the surface tension, 5i. Here we appeal to comparison of spontaneous 
decay to stimulated decay. The former dictates A'^ -C 1 in Eq.(15): the latter requires 
further discussion. 

Spontaneous dark energy decay brings us to the question of whether such decay can 
be initiated in an environment existing within our Universe. The question is analogous to 
one of electroweak phase transition in high energy particle collision. This was first raised 
in [30] and revisited for cosmic-ray collisions in [33]. That was in the context of the 
standard-model Higgs vacuum and the conclusion is that high-energy colliders are safe 
at all present and planned foreseeable energies because much more severe conditions 
have already occurred (without disaster) in cosmic -ray collisions within our galaxy. 
More recently, this issue has been addressed in connection with fears that the Relativistic 
Heavy Ion Collider (RHIC) might initiate a diastrous transition but according to careful 
analysis[34, 35] there was no such danger. 

The energy density involved for dark energy is some 58 orders of magnitude smaller 
[(lO^^eV)* compared to {300GeV)^~\ than for the electroweak case and so the nucle- 
ation scales are completely different. One is here led away from microscopic towards 
astronomical size scales. 

The energy density of Eq.(ll) is so readily exceeded that the critical radius cannot 
be microsopic. Think first of a macroscopic scale e.g. 1 meter and consider a magnetic 
field practically-attainable in bulk on Earth such as 10 Tesla. Its energy density is given 
by 

1 1 

Pmag - (16) 



Using the value jio = An x lO^'^NA^^ and IT = 6.2 x 10^^ (MeV.s.m^^) leads to an 
energy density p^ag — 2.5 x 10^'^ eV / {mm)^ , over 20 orders above the value of Eq.(ll) 
for the interior of nucleation. Magnetic fields in bulk exist in galaxies with strength 
~ ljuG and the rescaling by then would give p^ag ~ (2.5 x 10~^eV){mm)~^, slightly 
below the value in Eq.(l 1). 

Assuming the dark energy can exchange energy with magnetic energy density the ob- 
served absence of stimulated decay would then imply a critical radius of at least galactic 
size, say, ~ lOkpc. Using Eq(14) then gives for the surface tension Si > 10^^ (mm)~^ 
and number of nucleations in Eq.(15) < exp{—lO'^^). The spontaneous decay is thus 
glacial. Note that the dark energy has appeared only recently in cosmological time and 
has never interacted with background radiation of comparable energy density. Also, this 
nucleation argument does not require w < — 1. 



DISCUSSION 



As a first remark, since the critical radius for nucleation is astronomical, it appears 
that the instability cannot be triggered by any microscopic process. While it may be 
comforting to know that the dark energy is not such a doomsday phenomenon, it also 
implies at the same time the dreadful conclusion that dark energy may have no micro- 
scopic effect. If any such microscopic effect in a terrestrial experiment could be found, it 
would be crucial in investigating the dark energy phenomenon. We note that the present 
arguments are less model-dependent than those in [26]. 

In closing one may speculate how such stability arguments may evolve. One may 
expect most conservatively that the value w = — 1 will eventually be established empir- 
ically in which case both quintessence and the "phantom menace" will be irrelevant. In 
that case, indeed for any w, we may still hope that dark energy will provide the first con- 
nection between string theory and the real world as in e.g. [23]. Even if precise data do 
establish w < —1, as in the "phantom menace" scenario, the dark energy stability issue 
is still under control. 
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